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Abstract 

Using results on soft-collinear factorization for inclusive i?-meson decay distributions, a 
systematic study of the partial B — > decay rate with a cut > Eq on photon energy 
is performed. For values of Eq < 1.9 GeV the rate can be calculated without reference to 
shape functions. The result depends on three large scales: rrib, ^ym^,A, and A = mf,—2EQ. 
The sensitivity to the scale A 1.1 GeV (for Eq 1.8 GeV) introduces significant 
uncertainties, which have been ignored in the past. Our new prediction for the B Xg^ 
branching ratio with E^ > 1.8 GeV is Bi{B X^'y) = (3.44 ± 0.53 ± 0.35) x 10"^, 
where the errors refer to perturbative and parameter uncertainties, respectively. The 
implications of larger theory uncertainties for New Physics searches are explored with 
the example of the type-II two-Higgs-doublet model. 
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1 Introduction 



Given the prominent role of -B — X^'-) decay in searching for physics beyond the Standard 
Model, it is of great importance to have a precise prediction for its inclusive rate and CP 
asymmetry in the Standard Model. The total inclusive B Xs7 decay rate can be calculated 
using a conventional operator-product expansion (OPE) based on an expansion in logarithms 
and inverse powers of the 6-quark mass. However, in practice experiments can only measure 
the high-energy part of the photon spectrum, > Eq, where typically Eq = 2GeV or 
slightly below (measured in the 5-meson rest frame) ^1 12] . With E^ restricted to be close 
to the kinematic endpoint at Mb/2, the hadronic final state Xg is constrained to have large 
energy Ex ^ Mb but only moderate invariant mass Mx ~ (M^Aqcd)^''^- In this kinematic 
region, an infinite number of leading-twist terms in the OPE need to be resummed into a 
non-perturbative shape function, which describes the momentum distribution of the 6-quark 
inside the B meson H] . 

Conventional wisdom based on phenomenological studies of shape-function effects says 
these effects are important near the endpoint of the photon spectrum, but they can be ignored 
as soon as the cutoff Eq is lowered below about 1.9 GeV. In other words, there should be an 
instantaneous transition from the "shape-function region" of large non-perturbative correc- 
tions to the "OPE region", in which hadronic corrections to the rate are suppressed by at 
least two powers of AqcD/'^b- Below, we argue that this notion is based on a misconception. 
While it is correct that once the cutoff Eq is chosen below 1.9 GeV the decay rate can be 
calculated using a local short-distance expansion, we show that this expansion involves three 
"large" scales. In addition to the hard scale mj,, an intermediate scale i/mbA corresponding 
to the typical invariant mass of the hadronic final state X<j, and a low scale A = mi, — 2Eq 
related to the width of the energy window over which the measurement is performed, become 
of crucial importance. The precision of the theoretical calculations is ultimately determined 
by the value of the lowest short-distance scale A, which in practice is of order 1 GeV or only 
slightly larger. The theoretical accuracy that can be reached is therefore not as good as in 
the case of a conventional heavy-quark expansion applied to the B system. More likely, it 
is similar to (if not worse than) the accuracy reached, say, in the description of the inclusive 
hadronic decay rate of the r lepton. 

While we are aware that this conclusion may come as a surprise to many practitioners in 
the field of flavor physics, we believe that it is an unavoidable consequence of our analysis. 
Not surprisingly, then, we find that the error estimates for the B — >■ branching ratio 
that can be found in the literature are, without exception, too optimistic. Since there are 
unknown a^(A) corrections at the low scale A ~ 1 GeV, we estimate the present perturbative 
uncertainty in the B branching ratio with Eq in the range between 1.6 and 1.8 GeV 

to be of order 10-15%. In addition, there are uncertainties due to other sources, such as the 
h- and c-quark masses. The combined theoretical uncertainty is of order 15-20%, about twice 
as large as what has been claimed in the past. While this is a rather pessimistic conclusion, 
we stress that the uncertainty is limited by unknown, higher-order perturbative terms, not by 
non-perturbative effects, which we find to be under good control. Therefore, there is room for 
a reduction of the error by means of well-controlled perturbative calculations. 
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2 QCD factorization theorem 



Using recent results on the factorization of inclusive 5- meson decay distributions 0121, it is 
possible to derive a QCD factorization formula for the integrated B X^'y decay rate with a 
cut > Eq on photon energy. In the region of large Eq, the leading contribution to the rate 
can be factorized in the form [Zj 

^f'TjEo) = ^ \Vt,V;fml{f^,) \HM\' Urif,,, f..) (1) 

X / dP+{MB- P+Y dujmbJ{mb{P+-uj),Hi)S{iu,Hi), 
Jo Jo 

where A^; = Mb — 2Eq is twice the width of the window in photon energy over which the 
measurement of the decay rate is performed. The variable P+ = Ex — \Px\ is the "plus 
component" of the 4- momentum of the hadronic final state Xg, which is related to the photon 
energy by P+ = Mb — 2E^. The endpoint region of the photon spectrum is defined by the 
requirement that P+ < A^; ^ Mb, in which case is called a hard-collinear momentum jH]. 

In the factorization formula, /i/t ~ mi, is a hard scale, while /ij ~ y^m^jAqcD is an interme- 
diate hard-collinear scale of order the invariant mass of the hadronic final state. The precise 
values of these matching scales are irrelevant, since the rate is formally independent of fih and 
/ij. The hard corrections captured by the function H^{fih) result from the matching of the ef- 
fective weak Hamiltonian of the Standard Model (or any of its extensions) onto a leading-order 
current operator of soft-collinear effective theory (SCET) 0. At tree level, H^{fih) = Cj^{nh) 
is equal to the "effective" coefficient C^^ = Cj^ — | C5 — Cq. The expression valid at next- 
to- leading order can be found in Ref. jTj. The function H^{fih) is multiplied by the running 
6-quark mass rribifJ'h) defined in the MS scheme, which is part of the electromagnetic dipole 
operator Q^^. 

The jet function J{mb{P^ — uj),fii) in (jT)) describes the physics of the final-state hadronic 
jet. An expression for this function valid at next-to-leading order in perturbation theory has 
been derived in Refs. 13111 • The perturbative expansion of the jet function can be trusted as 
long as /if ~ nib A with A = rrib — 2Eq -C Mb- Note that the "natural" choices fih oc rrib and 
/if = rnfy fii with fti independent of remove all reference to the 6-quark mass (other than in 
the arguments of running coupling constants) from the factorization formula. 

The shape function S{u},fii) parameterizes our ignorance about the soft physics associated 
with bound-state effects inside the B meson 01^. Its naive interpretation is that of a parton 
distribution function, governing the distribution of the light-cone component of the residual 
momentum of the b quark inside the heavy meson. Once radiative corrections are included, 
however, a probabilistic interpretation of the shape function breaks down [S]. For convenience, 
the shape function is renormalized in at the intermediate hard-collinear scale /ij rather than 
at a hadronic scale /ihad- This removes any uncertainties related to the evolution from /ij to 
/ihad- Since the shape function is universal, all that matters is that it is renormalized at the 
same scale when comparing different processes. 

The last ingredient in the factorization formula is the function Ui{^h, ^i), which describes 
the renormalization-group (RG) evolution of the hard function l-ff-yp from the high match- 
ing scale fih down to the intermediate scale /ij, at which the jet and shape functions are 
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Figure 1: Dependence of the three scales fit = (sohd), /i, = y/nibA (dashed), and /xq = A 
(dash-dotted) on the cutoff Eq, assuming = 4.7 GeV. The gray area at the bottom shows 
the domain of non-perturbative physics. The hght gray band in the center indicates the region 
where the MSOPE must be apphed. 

renormahzed. The exact expression for this quantity and its perturbative expansion vahd at 
next-to-next-to- leading logarithmic order can be found in Ref. [7,. 

As written in (P), the decay rate is sensitive to non-perturbative hadronic physics via 
its dependence on the shape function. This sensitivity is unavoidable as long as the scale 
A = mb — 2EQ is a hadronic scale, corresponding to the endpoint region of the photon spectrum 
above, say, 2 GeV. Here we are interested in a situation where Eq is lowered out of the shape- 
function region, such that A can be considered large compared with Aqcd- For orientation, 
we note that with ruh = 4.7 GeV and the cutoff Eq = 1.8 GeV employed in a recent analysis 
by the Belle Collaboration j2] one gets A = 1.1 GeV. The values of the three relevant physical 
scales as functions of the photon-energy cutoff Eq are shown in Figure ^ This plot illustrates 
the fact that the transition from the shape-function region to the region where a conventional 
OPE can be applied is not abrupt but proceeds via an intermediate region, in which a short- 
distance analysis based on a multi-stage OPE (MSOPE) can be performed. The transition 
from the shape-function region into the MSOPE region occurs when the scale A becomes 
numerically (but not parametrically) large compared with Aqcd- Then terms of order a" (A) 
and (Aqcd/ A)"', which are non-perturbative in the shape-function region, gradually become 
decent expansion parameters. Only for very low values of the cutoff {Eq < 1 GeV or so) it is 
justified to treat A and \/mhA as scales of order m^. 

Separating the contributions associated with these scales requires a sophisticated multi- 
step procedure. The first step, the separation of the hard scale from the intermediate scale, 
has already been achieved in (|T)). To proceed further, we use that integrals of smooth weight 
functions with the shape function S{u,fi) can be expanded in a series of forward i?-meson 
matrix elements of local operators in heavy-quark effective theory (HQET) JH! , provided that 
the integration domain is large compared with Aqcd |S1 IE] • The perturbative expansions of 
the associated Wilson coefficient functions can be trusted as long as /i ~ A. In order to 
complete the scale separation, it is therefore necessary to evolve the shape function in (0) 
from the intermediate scale fii ~ y/rrihA down to a scale /xq ~ A. This can be achieved using 
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the analytic solution to the integro-differential RG evolution equation for the shape function 
in momentum space O [TT] . 

As a final comment, we stress that the main purpose of performing the scale separation 
using the MSOPE is not that this allows us to resum Sudakov logarithms. Indeed, the "large 
logarithm" ln(mb/A) ^ 1.5 is only parametrically large, but not numerically. What is really 
important is to disentangle the physics at the low scale fiQ ~ A, which is "barely perturbative" , 
from the physics associated with higher scales, where a short-distance treatment is on much 
safer grounds. The MSOPE allows us to distinguish between the three coupling constants 
asirrib) ^ 0.22, as{V^^) ^ 0.29, and a,(A) ^ 0.44 (for A = 1.1 GeV), which are rather 
different despite the fact that there are no numerically large logarithms in the problem. Given 
the values of these couplings, we expect that scale separation between A and nib is as important 
as that between rrih and the weak scale Mw- 



3 Calculation of the shape-function integral 

The scale dependence of the renormalized shape function is governed by an integro-differential 
RG evolution equation, whose exact solution in momentum space can be found using a tech- 
nique developed in Ref. jTJ. The result takes the remarkably simple form 



du' 



T{rj) Jo iil{u: - Cj'Y-^ ■ 



(2) 



The exact expression for the evolution function U2{^i, fio) can be found in Ref. [Zj, and 
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is given in terms of the cusp anomalous dimension [T^ . 

Relation Q accomplishes the evolution of the shape function from the intermediate scale 
down to the low scale /iq ~ A. The remaining task is to expand the integral over the shape 
function in (pj in a series of forward 5-meson matrix elements of local HQET operators 
of increasing dimension, multiplied by perturbative coefficient functions. This can be done 
whenever A = A^; — A = mb — 2Eq is large compared with Aqcd ,5j . The result is 



pleading ( rp \ 



32^ 



-lEV 



A 



X 



yjr^-l) (-Ai) 
2 3A2 



An 
+ ... 



4tt 



(4) 



where 



J{A) = 21n' 
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and h{ri) = ?/'(! + r]) +'~fE is the harmonic function generahzed to non-integer argument. Even 
though it is parametrically larger than ordinary power corrections of order [Aqcb / ^b)"^ , the 
"enhanced" Ai/A^ correction in (ji} remains small in the region of "perturbative" values of A, 
where the MSOPE can be trusted. The net effect amounts to a reduction of the decay rate 
by less than 5%. 

The rate in Q is formally independent of the three matching scales, at which we switch 
from QCD to SCET (/i/i), from SCET to HQET (/ij), and finally at which the shape-function 
integral is expanded in a series of local operators (fio). In practice, a residual scale dependence 
arises because we have truncated the perturbative expansion. Varying the three matching 
scales about their default values provides some information about unknown higher-order terms. 
In the limit where the intermediate and low matching scales fii and /iq are set equal to the hard 
matching scale fih, our result reduces to the conventional formula used in previous analyses of 
the B — > Xsj decay rate. However, this choice cannot be justified on physical grounds. 

In we have accomplished a complete resummation of (parametrically) large logarithms 
at next-to-next-to-leading logarithmic order in RG-improved perturbation theory, which is 
necessary in order to calculate the decay rate with 0{as) accuracy. Specifically, it means that 
terms of the form a'^L^ with k = {n — 1), . . . ,2n and L = ln(mfe/A) are correctly resummed 
to all orders in perturbation theory. To the best of our knowledge, a complete resummation 
at next-to-next-to-leading order has never been achieved before. Finally, we stress that the 
various next-to-leading order terms in the expression for the decay rate should be consistently 
expanded to order before applying our results to phenomenology. 

Up to this point, the 6-quark mass rrih entering the formula Q for the decay rate is 
defined in the on-shell scheme. While this is most convenient for performing calculations 
using heavy-quark expansions, it is well known that HQET parameters defined in the pole 
scheme suffer from infra-red renormalon ambiguities. It is necessary to replace them in favor 
of some physical, short-distance parameters. For our purposes, the "shape-function scheme" 
provides for a particularly suitable definition of the heavy-quark mass 0. The idea is that 
a good estimate of a shape-function integral can be obtained using the mean-value theorem, 
replacing u with 



Here m;,(A,/io) is the running shape-function mass, which depends on a hard cutoff A in 
addition to the renormalization scale /iq. The quantity A in the shape-function scheme 
is defined by the implicit equation A = A^; — A(A,/io) = mf,(A,/io) — 2Eq. The shape- 
function scheme provides a physical definition of mi,, which can be related to any other short- 
distance definition using perturbation theory. Based on various sources of phenomenological 
information including T spectroscopy and moments of inclusive S-meson decay spectra, the 
value of the shape-function mass at a reference scale /i* = 1.5 GeV has been determined as 




/c 




rrihifx,, IX,) = (4.65 ± 0.07) GeV 0. 
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The results discussed so far provide a complete description of the B Xg'j decay rate at 
leading order in the l/m^ expansion, where the two-step matching QCD ^ SCET — * HQET 
is well understood. For practical applications, however, it is necessary to include corrections 
arising at higher orders in the heavy-quark expansion. Most important are "kinematic" power 
corrections of order (A/mb)"", which are not associated with new hadronic parameters. Unlike 
the non-perturbative corrections, these effects appear already at first order in A/rrib, and they 
are numerically dominant in the region where A ^ Aqcd- Technically, the kinematic power 
corrections correspond to subleading jet functions arising in the matching of QCD onto higher- 
dimensional SCET operators, as well as subleading shape functions arising in the matching of 
SCET onto HQET operators. The corresponding terms are known in fixed-order perturbation 
theory, without scale separation and RG resummation [T8t I14j. To perform a complete RG 
analysis of even the first-order terms in A/rrib is beyond the scope of our discussion. Since for 
typical values of Eq the power corrections only account for about 15% of the B — >• Xgj decay 
rate, an approximate treatment suffices at the present level of precision. Details of how these 
corrections are implemented can be found in Ref . [7] . 

4 Ratios of decay rates 

The contributions from the three different short- distance scales entering the central result (j3)) 
and the associated theoretical uncertainties can be disentangled by taking ratios of decay rates. 
Some ratios probe truly short- distance physics (i.e., physics above the scale fih ~ rnt) and so 
remain unaffected by the new theoretical results presented above. For some other ratios, the 
short- distance physics associated with the hard scale cancels to a large extent, so that one 
probes physics at the intermediate and low scales, irrespective of the short-distance structure 
of the theory. 

Ratios insensitive to low-scale physics: Physics beyond the Standard Model may affect 
the theoretical results for the B — > Xs'j branching ratio and CP asymmetry only via the 
Wilson coefficients of the various operators in the effective weak Hamiltonian. As a result, 
the ratio of the B Xg'y decay rate in a New-Physics model relative to that in the Standard 
Model remains largely unaffected by the resummation effects studied in the present work. 
From Q, we obtain 



The power corrections would introduce some mild dependence on the intermediate and low 
scales i^i and /io, as well as on the cutoff Eq. 

Another important example is the direct CP asymmetry in B —>■ Xgj decays, for which we 
obtain 
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where is obtained from by CP conjugation, which in the Standard Model amounts to 
replacing the CKM matrix elements by their complex conjugates. It follows that predictions 
for the CP asymmetry in the Standard Model and various New Physics scenarios ^H] remain 
largely unaffected by our considerations. 



Ratios sensitive to low-scale physics: The multi-scale effects studied in this work result 
from the fact that in practice the B —>■ Xs'y decay rate is measured with a restrictive cut on 
the photon energy. These complications would be absent if it were possible to measure the 
fully inclusive rate. It is convenient to define a function F{Eq) as the ratio of the B — > Xsj 
decay rate with a cut Eq divided by the total rate, 

F(E,) ^ I'-'-fjl . (9) 



Because of a logarithmic soft-photon divergence for very low energy, it is conventional 
to define the "total" inclusive rate as the rate with a very low cutoff = mb/20. The 
denominator in the expression for F{Ef)) can be evaluated using the standard OPE, which 
corresponds to setting all three matching scales equal to fih- The numerator is given by our 
expression in (0}, supplemented by power corrections. 

Another important example of a ratio that is largely insensitive to the hard matching 
contributions is the average photon energy (E^), which has been proposed as a good way to 
measure the 6-quark mass ^Hl- The impact of shape-function effects on the theoretical predic- 
tion for this quantity has been investigated and was found to be significant Here we 
study the average photon energy in the MSOPE region, where a model-independent prediction 
can be obtained. It is structurally different from the one found using the conventional OPE in 
the sense that contributions associated with different scales are disentangled from each other. 
We stress that the hard scale fih ~ affects the average photon energy only via second-order 
power corrections. This shows that it is not appropriate to compute the quantity (E^) using 
a simple heavy-quark expansion at the scale rrif,, which is however done in the conventional 
approach jTHj. This observation is important, because information about moments of the 
B Xs'j photon spectrum is sometimes used in global fits to determine the CKM matrix el- 
ement \Vcb\- Keeping only the leading power corrections, which is a very good approximation, 
we find that {E^) only depends on physics at the intermediate and low scales fii and fiQ. For 
Eo = 1.8 GeV, we obtain (E^) ^ [2.222 + 0.254:as{VmA) + 0.009a,(A)] GeV ^ 2.30 GeV. 



5 Numerical results 

We are now ready to present the phenomenological implications of our findings. A complete 
list of the relevant input parameters and their uncertainties is given in Ref . [7] , where we also 
explain our strategy for estimating the perturbative uncertainty as well as the uncertainty due 
to parameter variations. 

We begin by presenting predictions for the CP-averaged B Xs'j branching fraction with 
a cutoff E^ > Eq applied on the photon energy measured in the i?-meson rest frame. Lowering 
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Eq below 2 GeV is challenging experimentally. The first measurement with = 1.8 GeV has 
recently been reported by the Belle Collaboration f7. It yields^ 



Br(5^X,7) = (3.38 ±0.30 ±0.29) ■ 10 

-E-/>1.8 GeV 
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(E^) = (2.292 ± 0.026 ± 0.034) GeV . (10) 

^ ^' Eo=1.8GeY ^ ' ^ ' 

For Eq = 1.8 GeV we have A ^ 1.1 GeV, which is sufficiently large to apply the formalism 
developed in the present work. (For comparison, the value Eq = 2.0 GeV adopted in the CLEO 
analysis implies A ^ 0.7 GeV, which we believe is too low for a short-distance treatment.) 
We find 

Br(5 ^X,7)|^ -isGeV " (3-44 ±0.53 [pert.] ±0.35 [pars.]) x 10'^, (11) 

where the ffist error refers to the perturbative uncertainty and the second one to parameter 
variations. The largest parameter uncertainties are due to the b- and c-quark masses. Our 
result is in excellent agreement with the experimental value shown in (|10p . Comparing the 
two results, and naively assuming Gaussian errors, we conclude that 

Br(E^X,7)exp-Br(5^X,7)sM<1.4-10-^ (95% CL) . (12) 

Mainly as a result of the enlarged theoretical uncertainty, this bound is much weaker than the 
one derived in Ref. where this difference was found to be less than 0.5 ■ 10~^. Hence, we 
obtain a much weaker constraint on New Physics parameters. For instance, for the case of the 
type-II two-Higgs-doublet model, we may use the analysis of Ref. [11] to deduce 



mH+ > (slightly below) 200 GeV (95% CL) , (13) 

which is significantly weaker than the constraints mH+ > 500 GeV (at 95% CL) and mH+ > 
350 GeV (at 99% CL) found in Ref. 

The function F{Eq) provides us with an alternative way to discuss the effects of imposing 
the cutoff on the photon energy. In contrast to the branching ratio, it is independent of 
several input parameters (e.g., mb(mfe), |Vj*Vtb|, tb, ^1,2), and it shows a very weak sensitivity 
to variations of the remaining parameters. We obtain 

F(1.8GeV) = (92+1^ [pert.] ± 1 [pars.])%. (14) 

This is the ffist time that this fraction has been computed in a model independent way. The 
result may be compared with the values (95.8 12^)% and (95 ± 1)% obtained from two studies 
of shape-function models ITT] , in which perturbative uncertainties have been ignored. We 
obtain a significantly smaller central value with a much larger uncertainty. 

The last quantity we wish to explore is the average photon energy. As mentioned above, 
this observable is very sensitive to the interplay of physics at the intermediate and low scales. 
The study of uncertainties due to parameter variations exhibits that the prime sensitivity is 



^To obtain the first result we had to undo a theoretical correction accounting for the effects of the cut 
> 1.8 GeV, which had been applied to the experimental data. 
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to the 6-quark mass, which is expected, since {E^) = m{,/2 + ... to leading order. The next- 
important contribution to the error comes from the HQET parameter Ai. To a very good 
approximation, we have 

{E,) =(2.27+°o°?)GeV+^-|^, (15) 

^ ^' £0=1.8 Gev ^ u-'J^^ 2 Arub 

where the error accounts for the perturbative uncertainty. The quantities 6mb and 6X1 pa- 
rameterize possible deviations of the relevant input parameters from their central values 
rrih = 4.65 GeV and Ai = — 0.25GeV^. Our prediction is in excellent agreement with the 
Belle result in flTUj) . This finding provides support to the value of the 6-quark mass in the 
shape-function scheme extracted in Ref. 0. We stress, however, that the large perturbative 
uncertainties in the formula for (E^) impose significant limitations on the precision with which 
mft can be extracted from a measurement of the average photon energy. Our estimate above 
implies a perturbative uncertainty of 5m;,[pert.] = Ij^QgMeV. This is in addition to twice the 
experimental error in the measurement of (E^), which at present yields 5mfe[exp.] = 86MeV. 



6 Conclusions and outlook 

We have performed the first systematic analysis of the inclusive decay B — > Xgj in the presence 
of a photon-energy cut E^ > Eq, where Eq is such that A = mi, — 2Eq can be considered large 
compared to Aqcd, while still A <^ m^. This is the region of interest to experiments at the B 
factories. The first condition (A ^ Aqcd) ensures that a theoretical treatment without shape 
functions can be applied. However, the second condition (A -C rrib) means that this treatment 
is not a conventional heavy-quark expansion in powers of as{mb) and Aqcd/'^t.;,. Instead, we 
have shown that three distinct short- distance scales are relevant in this region. They are the 
hard scale nib, the hard-coUinear (or jet) scale ^/mbA, and the low scale A. To separate the 
contributions associated with these scales requires a multi-scale operator product expansion 
(MSOPE). 

Our approach allows us to study analytically the transition from the shape-function re- 
gion, where A ~ Aqcd, into the MSOPE region, where Aqcd ^ A ^ m^, into the region 
A = 0{mb), where a conventional heavy-quark expansion applies. This is a significant im- 
provement over previous work. For instance, it has sometimes been argued that exactly where 
the transition to a conventional heavy-quark expansion occurs is an empirical question, which 
cannot be answered theoretically. Our formalism provides a precise, quantitative answer to this 
question. In particular, for B X^'y with a realistic cut on the photon energy, one is not in a 
situation where a short- distance expansion at the scale mb can be justified. The analysis makes 
it evident that the precision that can be achieved in the prediction of the B Xg'-f branching 
ratio is, ultimately, determined by how well perturbative and non-perturbative corrections can 
be controlled at the lowest relevant scale A, which in practice is of order 1 GeV. Consequently, 
we estimate much larger theoretical uncertainties than previous authors. These uncertainties 
are dominated by yet unknown higher-order perturbative effects. Non-perturbative, hadronic 
effects at the scale A appear to be small and under control. 

This is not the first time in the history of i? — > Xg'j calculations that issues of scale setting 
have changed the prediction and error estimate for the branching ratio (see, e.g., the discussion 
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in Ref. ^H|)- In our case, however, the change in perspective about the theory oi B Xg'j 
decay is more profound, as it imposes hmitations on the very vahdity of a short- distance 
treatment. If the short-distance expansion at the scale A fails, then the rate cannot be calcu- 
lated without resource to non-perturbative shape functions, which introduces an irreducible 
amount of model dependence. In practice, while A ^ 1.1 GeV (for Eq ^ 1.8 GeV) is probably 
sufficiently large to trust a short-distance analysis, it would be unreasonable to expect that 
yet unknown higher-order effects should be less important than in the case of other low-scale 
applications of QCD. 

Obtaining a precise prediction for the B X^'f decay rate in the Standard Model is an 
important target of heavy-flavor theory. The present work shows that the ongoing effort to 
calculate the dominant parts of the next-to-next-to-leading corrections in the conventional 
heavy-quark expansion is only part of what is needed to achieve this goal. Equally impor- 
tant will be to compute the dominant higher-order corrections proportional to a^(A) and 
a^(\/mbA), and to perform a renormalization-group analysis of the leading kinematic power 
corrections of order A /nib. In fact, our error analysis suggests that these effects are potentially 
more important that the hard matching corrections at the scale nib. 
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